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ABSTRACT 


The  purpose  of  this  thesis  is  to  develop  a  means  by  which 
the  transient  response  of  a  servomechanism  may  be  rapidly 
evaluated.  The  approach  consists  of  developing  a  method  of 
rapidly  determining  the  value  of  several  transient  response 
parameters  so  that  the  system  transient  response  may  be  sketched 
The  parameters  selected  for  this  purpose  are  time  and  magnitude 
of  peak  overshoot,  rise  time  and  settling  time. 

To  permit  evaluation  of  these  parameters  systems  are 
classified  according  to  the  number  and  type  of  singularities 
in  their  closed  loop  performance  functions.   Systems  considered 
are  underdamped  systems  with  two  or  three  singularities ,   For 
each  class  of  system,  loci  of  constant  values  of  overshoot  and 
constant  values  of  the  three  time  parameters  expressed  in  terms 
of  an  artificial  time  are  mapped  onto  a  complex  plane.  Para- 
meter loci  for  systems  with  only  two  singularities  are  mapped 
directly  onto  the  s -plane.   Loci  for  systems  with  three  singu- 
larities are  mapped  onto  an  auxiliary  complex  plane. 

The  values  of  the  four  transient  response  parameters  for 
any  system  of  a  particular  class  can  be  read  from  the  families 
of  loci  for  that  class.  The  magnitude  of  overshoot  can  be 
used  directly  and  the  three  time  parameters  can  be  easily  eon- 
verted  into  real  time.  The  initial  slope  of  the  transient 
response  and  the  slope  at  the  time  of  peak  overshoot  are  either 
known  or  easily  determined.   Also  the  transient  oscillating 
frequency  is  obtainable  directly  from  the  system  closed  loop 
performance  function,,   Using  this  information  the  transient 
response  of  a  particular  system  can  be  sketched  with  considerable 
accuracy.  This  sketch  can  then  be  used  as  a  basis  for  rapid 
system  evaluation. 

Thesis  Supervisor:   Robert  K.  Mueller 

Title:   Associate  Professor  of  Aeronautics  and  Astronautics 
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Chapter  I 
INTRODUCTION 

In  the  past  two  decades  the  importance  of  servome onanisms 
and  automatic  control  systems  in  general  has  increased  tre- 
mendously. Accordingly,  much  effort  has  been  expended  in  de- 
veloping methods  of  analysis  and  design.  The  result  is  an 
abundance  of  ideas  and  techniques  available  to  the  designer 
today.  Most  of  these  can  be  found  in  any  one  of  a  number  of 
current  automatic  control  system  textbooks. 

Basically  there  are  two  approaches  to  analysis  and  design, 
the  frequency  response  approach  and  the  time  domain  approach. 
Either  of  these  techniques  may  be  used  with  some  facility  in 
designing  systems  to  meet  both  time  domain  and  frequency  re- 
sponse specifications.   However,  if  time  domain  specifications 
are  to  be  met,  a  time  domain  analysis  must  eventually  be  made. 
This  is  true  regardless  of  which  technique  is  used  in  the 
actual  design. 

The  usual  time  domain  specification  is  one  involving 
transient  response,  that  is,  the  time  response  of  the  system 
to  a  suddenly  applied  input  signal.  The  most  common  input 
signal  used  is  the  unit  step.  Hereafter  we  will  use  the  term 
transient  response  exclusively  to  mean  the  time  response  of  a 
system  to  a  unit  step  input  signal  with  all  initial  conditions 
equal  to  zero. 
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The  problem  of  obtaining  the  transient  response  of  servo- 
mechanisms  is  one  of  considerable  magnitude.  The  transient 
response  of  linear  second  order  systems  has  been  studied  ex- 
tensively and  is  well  known.  The  transient  responses  of  high- 
er order  systems  and  second  order  systems  that  have  been 
complicated  by  addition  of  a  zero  are  more  difficult  to  obtain 
and  are  less  well  known.  The  exact  transient  response  of  any 
linear  system  can  always  be  obtained  by  factoring  the  charac- 
teristic equation,  solving  for  the  output  as  a  function  of 
time  and  plotting  this  function.  However,  this  procedure  in- 
volves considerable  labor  and  is  quite  time  consuming. 

One  method  of  attempting  to  avoid  such  a  calculation  is  to 
consider  the  complex  roots  of  the  characteristic  equation  to 
be  dominant  and  use  a  second  order  approximation  of  the  more 
complicated  system.  As  the  value  of  the  real  roots  and  zeros 
become  smaller  with  respect  to  the  undamped  natural  frequency 
of  the  system  this  approximation  ceases  to  be  valid.  When  the 
second  order  approximation  ceases  to  be  valid  another  approxi- 
mation developed  by  Chu  (l)  may  be  employed.  This  method  re- 
duces the  labor  involved  in  obtaining  the  transient  response 
by  applying  corrections  to  the  second  order  equations  for 
zeros  and  other  roots  of  a  system.   Unfortunately  this  approxi- 
mation also  becomes  invalid  as  the  values  of  the  real  roots 
and  zeros  become  even  smaller. 

If  the  frequency  response  approach  is  used,  obtaining  the 
transient  response  is  usually  more  difficult  because  the  system 
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characterise  Ration  is  „ot  known  ln  factored  fom_ 

b  i7rphioai  intesration  schemes  such  - -  *°— 

tL  ^  S°meUmeS  US6d  t0  °btal"  -  <™»>ate 

talent  response.     These  methods  also  ^^^ 

labor  and  time.  «aDie 

In  this  thesis  we  restrict  +**>  ** 

StriCt  tne  discussion  to  linear 

I  r1"  eXhibU  ^  Pr°Perty  ■"-**•*■*.  -at  u. 

a :  t  seor  rch  the  steady  state — *--  -* — 

staT lli         3l6nal  "  aPPll0d-  NOtl°e  «*  «-»*• 
stability  is  Implied. 

With  this  restriction  „e  can  mafce  some  general  obser- 
vations about  tne  transient  response  of  a„y  such  system. 

a-  At  zero  time  the  value  of  the  output 
is  zero. 

°.  At  a  sufficiently  large  time  the 
value  of  the  output  becomes  unity. 

o.  For  critically  damped  and  over- 
damped  systems  the  output  increases 
monotonically  with  time  until  it 
reaches  its  steady  state  value  of 
unity.  This  case  is  of  little 
practical  Interest. 

«•  A  more  common  and  more  practical 
case  is  that  of  an  underdamped 
astern,  m  this  case  the  output 
rises  to  a  value  greater  than 
unity,  and  may  or  may  not  oscillate 
about  unity  before  assuming  its 
steady  state  value. 
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Based  on  the  above  assumptions  and  observations  we  may 
define  the  following  transient  response  parameters: 

a.  Peak  Overshoot  -  M       The  amount  by  which  the  maxi- 
mum value  of  the  output  exceeds  the  steady  state  value. 

b.  Time  of  Peak  Overshoot  -  t     The  time  at  which  the 

E 

maximum  value  of  output  occurs. 

c.  Rise  Time  -  t     The  time  at  which  the  value  of  the  out- 
put is  unity  for  the  first  time. 

d.  Settling  Time  -  t     The  mathematical  expression  for  the 
transient  response  of  systems  which  have  only  one  pair 

of  complex  roots  may  consist  of  three  parts.  These  are 
a  constant  equal  to  the  steady  state  value  of  output, 
decaying  exponential  terms,  and  a  damped  sinusoidal 
term.  The  settling  time  is  the  time  at  which  the  sum 
of  the  error  due  to  the  exponential  terms  plus  the 
height  of  the  envelope  of  the  damped  sinusoid  is  equal 
to  two  per  cent  of  the  steady  state  output  value. 

A  sketch  of  a  typical  transient  response  with  the  above- 
defined  parameters  indicated  is  shown  in  Figure  1-1. 

If  the  four  transient  response  parameters  defined  above 
and  the  transient  oscillating  frequency  of  a  system  are  known, 
a  sketch  which  very  closely  approximates  the  actual  transient 
response  can  be  constructed.  A  means  of  rapidly  determining 
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TYPICAL  TRANSIENT  RESPONSE 
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these  parameters  would  then  be  a  valuable  aid  to  the  de- 
signer in  evaluating  the  transient  performance  of  a  system 
without  computing  the  exact  transient  response. 

In  this  thesis  a  technique  which  permits  rapid  graphical 
determination  of  these  parameters  by  inspection  is  developed 
for  certain  classes  of  systems.  The  technique  involves  the 
mapping  of  lines  of  constant  value  of  these  parameters  onto 
a  complex  plane  for  a  particular  class  of  system.  The  values 
of  these  parameters  for  all  systems  of  that  class  can  then 
be  read  directly  from  the  complex  plane  plots. 
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Chapter  II 
BASIS  OF  THE  MAPPING  TECHNIQUE 

The  determination  of  the  transient  response  parameters  of 
a  given  system  is  most  easily  approached  through  the  per- 
formance function,  the  singularities  of  which  can  be  displayed 
on  a  complex  plane.  In  this  approach  the  Laplace  transform 
with  the  associated  complex  function  theory  is  used  as  a 
mathematical  tool. 

The  mapping  technique  developed  in  thds  thesis  is  based 
on  the  closed  loop  performance  function,  the  poles  of  which 
are  actually  the  roots  of  the  system  characteristic  equation. 
The  technique  as  developed  is  applicable  to  systems  whose 
closed  loop  performance  functions  contain  one  and  only  one 
pair  of  complex  conjugate  poles. 

Now  consider  such  a  performance  function  consisting  of 
one  pair  of  complex  poles  and  an  assortment  of  real  poles 
and  zeros.  The  poles  and  zeros  of  the  system  can  be  plotted 
on  a  complex  plane.  If  the  locations  of  the  real  poles  and 
zeros  are  fixed,  a  set  of  values  for  the  transient  response 
parameters  defined  in  Chapter  I  can  be  determined  for  each 
possible  position  of  the  complex  poles.  If  the  values  of 
each  parameter  are  determined  for  a  sufficient  number  of 
complex  pole  locations,  the  loci  of  constant  values  of  es 
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parameter  can  then  be  drawn  on  the  complex  plane.  These  loci 
are  functions  of  the  complex  pole  locations  but  are  valid 
only  for  the  fixed  configuration  of  real  poles  and  zeros. 

A  simple  example  of  the  loci  of  constant  values  of  a 
parameter  exists  in  the  case  of  a  pure  second  order  system, 
that  is,  a  system  for  which  the  closed  loop  performance 
function  consists  of  ohe  pair  of  complex  poles  and  no  real 
poles  or  zeros.  In  this  case  the  overshoot  is  a  function 
of  damping  ratio  only  and  the  lines  of  constant  damping  ratio 
are  also  lines  of  constant  values  of  overshoot.  These  are 
radial  lines  emanating  from  the  origin  of  the  complex  plane. 
This  relationship  is  often  displayed  as  a  plot  of  overshoot 
versus  damping  ratio  as  is  shown  by  Thaler  (5). 

Clearly  each  configuration  of  real  poles  and  zeros  re- 
sults in  a  different  set  of  loci.  For  example,  the  addition 
of  a  real  pole  or  zero  to  the  second  order  case  just  con- 
sidered results  in  lines  of  constant  overshoot  which  ar©  not 
straight  lines.  In  fact,  the  exact  location  and  shape  of 
the  lines  of  constant  overshoot  are  determined  by  the  actual 
locations  of  the  real  pole  or  zero. 

This  leads  to  the  need  for  a  method  of  classifying 
systems  according  to  the  real  pole  and  zero  configuration. 
In  order  to  facilitate  the  investigation  of  the  sets  of 
loci  associated  with  different  systems,  the  systems  have 
been  classified  as  follows: 
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a.   Pure  Second  Order  System  -  Closed  loop 
performance  function  consists  of  two 
complex  poles  only. 

b0  Second  Order  System  With  One  Zero  - 
Closed  loop  performance  function  con- 
sists of  one  pair  of  complex  poles  and 
one  real  zero. 

c.  Pure  Third  Order  System  -  Closed  loop 
performance  function  consists  of  one 
pair  of  complex  poles  and  one  real  pole. 

These  three  classes  of  systems  have  been  nondimen- 
sionalized  so  that  one  set  of  loci  of  constant  values  of 
parameters  is  sufficient  for  all  systems  of  that  particular 
class.  The  details  of  the  nondimensionalizing  process  and 
the  complex  plane  plots  of  loci  of  constant  parameter  values 
are  shown  in  Chapters  III,  IV,  and  V.  Use  of  the  families 
of  loci  determined  in  these  chapters  requires  that  the  closed 
loop  performance  function  of  a  system  be  known  in  factored 
form. 

Only  the  second  quadrant  of  the  complex  plane  is  shown 
for  two  reasons. 

a.  Only  stable  systems  are  considered  so 
the  poles  are  constrained  to  be  in  the 
second  and  third  quadrants. 
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b.  The  complex  poles  are  conjugates  so 
the  third  quadrant  loci  are  images  of 
the  second  quadrant  loci. 

In  Chapter  VI  the  technique  is  extended  to  third  order 
systems  with  one  or  two  zeros  by  means  of  an  algebraic 
manipulation  which  permits  evaluating  the  parameters  of 
two  or  more  systems  of  the  classes  listed  above  and  com- 
bining the  results. 
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Chapter  III 
THE  PURE  SECOND  ORDER  SYSTEM 

The  second  order  system  Is  the  simplest  type  to  which 
the  previously  outlined  mapping  technique  is  applicable. 
Its  characteristics  have  been  extensively  studied  and  almost 
every  servomechanism  text  contains  a  complete  description  of 
its  transient  response.  The  details  of  this  transient  re- 
sponse will  not  be  repeated  here. 

In  this  chapter  we  are  concerned  with  what  may  be  called 
a  pure  second  order  system,  that  is,  one  whose  closed  loop 
performance  function  consists  of  only  two  complex  conjugate 
poles .  The  performance  function  of  such  a  system  is 

pp(s)  -  HH  -  ;  ^  X  h2  2  ft) 

R(s;   (s+ar  +  *> 

The  poles  of  this  system  are  shown  on  the  s -plane  in 
Figure  III-l. 

When  R(t)  is  the  unit  step  function  the  transient  re- 
sponse of  the  system  becomes 


2    2 
C(t)  =  1  +   1  a \ +  b   e"at  sin(bt-O)  (2) 


^b" 

where  0=  arc  tan  — —■ 

— s 
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FIGURE  III-l 
POLES  OF  THE  PURE  SECOND  ORDER 
SYSTEM  ON  THE  8 -PLANE 
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This  is  a  very  familiar  expression,  although  it  is  usually 
written  in  polar  coordinates,  i.e.,  in  terms  of  damping 
ratio  and  undamped  natural  frequency. 

It  is  convenient  to  define  an  artificial  time  unit 
as 

T  -  |a|t  (3) 

The  expression  for  the  transient  response  in  terms  of  T 
becomes 


C  (T )  -  1  + 


2 
%  +  1  e"T  sin(fr  -  9)  (4) 


In  terms  of  the  artificial  time,  T,  the  entire 
transient  response  is  determined  if  the  ratio  of  a  to  b 
is  known.  This  is  equivalent  to  knowing  the  damping  ratio. 

Time  and  Magnitude  of  Peak  Overshoot 

The  slope  of  the  transient  response  is 

C(t)  =  a  b+b  e"at  sin  bt  (5) 


or 


.2  ,  v2 


0(T)  ~  a  b+  b  e"T  sin|r  (6) 


From  the  nature  of  the  transient  response  it  can  be 

seen  that  the  peak  overshoot  occurs  at  the  first  time  after 

bT 
time  zero  that  the  slope  is  zero.  This  requires  that  sin  ~- 

be  equal  to  zero,  or 
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|r  =  «  (7) 

Therefore    T  =  *#  (8) 


In  real  time 


t  -  J  -  5  (9) 

p   iai   b  v*' 


as  is  shown  by  Truxal  (4). 

The  loci  of  constant  values  of  time  of  peak  overshoot 
(T  )  are  shown  mapped  onto  the  s-plane  in  Figure  III-2. 
Although  a  scale  is  shown  on  this  plot  for  convenience,  the 
loci  are  independent  of  scale  as  they  depend  only  on  the 
ratio  of  a  to  b.   Hence  any  convenient  scale  may  be  sub- 
stituted for  the  scale  shown. 

The  magnitude  of  peak  overshoot  can  be  determined  by 
substituting  T  into  the  expression  for  C(T),  Equation  4. 
The  result  of  this  substitution  is 

C(T  )  =  1  +  e   b  (10) 

The  magnitude  of  peak  overshoot  is  the  value  of  the  ex- 
ponential term  and  once  more  dependent  only  on  the  ratio 

a 

F*  _  ra 

M  -  100  e  ""*>  %  (11) 

Loci  of  constant  magnitude  of  peak  overshoot  are  shown 
mapped  onto  the  s-plane  in  Figure  III -5.  They  display  the 
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II 


same  information  that  is  usually  presented  in  the  form  of 
overshoot  versus  damping  ratio  plots,  as  has  been 
mentioned. 

Rise  Time 

Rise  time  has  been  defined  as  the  first  time  after  zero 
time  when  the  value  of  the  transient  response  is  unity.  In- 
spection of  Equation  4  shows  that  this  occurs  when 


sin &T  -  9)  =*  0  (12) 

a 


or 


§T  -  0  -  0  (15) 

Therefore     Tr  "  f°  =  f  arctan  1£  (l20 


Loci  of  constant  values  of  artificial  rise  time  (T  ) 
are  shown  mapped  onto  the  s -plane  in  Figure  III-4. 


Settling  Time 

Settling  time  has  been  defined  as  the  time  when  the 
height  of  the  transient  envelope  becomes  two  percent  of 
the  steady  state  output  value.  This  occurs  when 


a2       -T 

±2  +  1  e  1  =  .02  (15) 
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Therefore 

t  =  -In  .02  +  i  ln(\  +  i)  (16) 

s  ^    b^ 

Loci  of  constant  values  of  artificial  settling  time 
(T_)  are  shown  mapped  onto  the  s -plane  in  Figure  III -5. 

5 


Use  of  the  Curves 

To  use  the  curves  shown  in  Figures  III -2  through  III -5, 
it  is  necessary  to  determine  a  and  b.  Effectively  this 
means  that  the  roots  of  the  system  characteristic  equation 
must  be  known.  The  method  used  in  obtaining  these  roots  is 
immaterial . 

Once  the  roots  are  known  some  convenient  arbitrary 

scale  can  be  chosen  for  the  loci  plots.  Since  these  are 

complex  plane  plots  the  real  and  imaginary  scales  must  be 

identical.  The  point  -a  +  jb  can  then  be  located  on  each 

loci  plot  and  the  values  of  M  .  T  .  T  ,  and  T„,  read  off 

p   p   r      s 

directly. 

The  three  artificial  time  parameters  can  then  be  con- 
verted to  real  times  by  dividing  each  by  the  value  of  a. 
The  transient  response  may  then  be  rapidly  sketched  by 
marking  off  known  points,  remembering  that: 

a.  The  slope  is  zero  at  t  *  0 
and  t  »  t  . 
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b.  The  transient  oscillating  frequency 
is  ■£-  cycles  per  second. 

c.  The  sine  function  is  zero  at  t  -  t 

r 

and  increasing  with  time. 

Figure  III -6  is  an  example  of  the  procedure  for  the 
performance  function  shown.  An  exact  curve  is  included 
for  comparison. 

Conclusion 

It  has  been  shown  that  the  loci  of  constant  values  of 
T  ,  T  ,  Ts,  and  M  can  be  mapped  directly  onto  the  s-plane 
for  the  pure  second  order  system.  In  addition,  these  four 
parameters  and  the  transient  oscillating  frequency  present 
enough  information  about  the  transient  response  of  such  a 
system  to  permit  a  reasonably  accurate  sketch  to  be  made. 
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Chapter  IV 
THE  SECOND  ORDER  SYSTEM  WITH  ONE  ZERO 

The  addition  of  a  real  zero  to  a  pure  second  order  system 
results  in  a  different  type  of  second  order  system  whose 
closed  loop  performance  function  is 

PF(S)  =  *2  +*>2  .  (g  +  d)  (1) 

d        (s+ar  +  *> 

The  singularities  of  this  function  are  shown  on  the  s-plane 
in  Figure  IV-1.  The  constant  which  is  usually  interpreted 
as  system  gain  reflects  the  assumption  of  unity  feedback 
characteristics . 

As  pointed  out  by  Truxal  (4),  the  Laplace  transform  of 
the  response  of  such  a  system  to  a  unit  step  input  may  be 
written  as 

c(s)4   a2+f   g     +i    *2y>2  g     (2) 

s       (s+a)2  +  b2  d       (s+a)*  +  bd 

This  response  is  the  sum  of  the  unit  step  response  of  a  pure 
second  order  system  plus  a  portion  of  the  impulse  response 
of  the  same  system.  The  addition  of  the  impulse  response 
makes  the  system  relatively  more  oscillatory  than  would  be 
the  case  without  the  zero.  The  magnitude  of  this  effect 
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FIGURE  IV  -1 

POLES  AND  ZERO  ON  THE  3 -PLANE 

Second  Order  System  with  One  Zero 
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depends  on  the  position  of  the  zero  relative  to  the 
position  of  the  complex  poles . 

It  is  also  useful  to  regard  the  addition  of  the  zero 
as  a  consequence  of  the  use  of  proportional  plus  derivative 
control  in  the  forward  path  of  a  conventional  unity  feed- 
back servomechanisra.  Another  equally  plausible  view  is  to 
consider  this  case  as  an  approximation  of  a  higher  order 
system,  that  is,  one  in  which  the  singularities  that  are 
far  out  in  the  left  half -plane  are  neglected. 

The  expression  for  the  transient  response  of  a  system 
whose  performance  function  is  given  by  Equation  1  is 


c(t) 


\  2   2   I 

1  +  ^a  +db   \|  (d-a)2  +  b2  e"at  sin(bt  +  0   -  0)    (3) 


where  0  -  arctan  a—-,  and  0  =  arctan  -~ 

r  d-a  -a 


The  slope  of  this  transient  response  is 


.2  .  .2 


,2   ,2  .-at 


C(t)  -  ^         \(d-ar  +  b  *        sln(bt  +  ft)  (4) 


We  may  now  express  a  and  b  in  terms  of  d  by  defining 

mi  =  itr  (5) 

m2  -  iSr  (6) 

Figure  IV-2  illustrates  these  quantities. 
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FIGURE  IV -2 


ILLUSTRATION  OF  THE  QUANTITIES  m1  AND  m2 
Second  Order  System  with  One  Zero 
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Substituting  these  relationships  into  the  expressions 
for  transient  response  and  slope  of  transient  response, 
we  have 


C(t) 


\  m   +  nu  \(l-m, )   +  nu   -m  dt 

1  + ± £ ± £_  e  x 

nu 


sin(m2dt-^-0)    (7) 


m. 


2     2 

C(t)  =  ^-A_ Z-L  ^(i-n^)2  +  m22  e  1     sin(m2dt-^) 


(8) 


where 


0   =  arctan 


9  -  arctan 


m. 


l-m„ 


m. 


-m. 


We  now  define  an  artificial  time  unit  for  this  type  of 
system  as 

T  =  Idl  t 


In  terms  of  the  artificial  time,  T,  the  expressions  for 
transient  response  and  slope  of  transient  response  become 


C(T) 


C(T) 


K2  +m22  ^(1-m^2  +  m22  -n^T 


sin(moT^-0)    (9) 


dTi 


mn   +  m0  ) 


m. 


~^~      1  2  -m  T 
(1-m  )  +  m2  e  >> '  3in  (m  T-$ ) 


(10) 


where  0  and  0   are  defined  as  before. 
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It  is  important  that  the  reader  visualize  what  has  been 
done.  The  simplest  view  is  that  the  transient  response  is 
expressed  as  a  function  of  an  arbitrary  time,  T,  and  that 
this  response  depends  only  on  the  position  of  the  complex 
poles  relative  to  the  zero.  One  merely  has  to  adopt  the 
magnitude  of  d  as  the  unit  of  measure  and  determine  m_  and 
m?  by  measuring  a  and  b  on  the  s -plane  using  this  new  unit. 
In  terms  of  T  the  transient  response  is  then  independent  of 
the  actual  numbers  associated  with  a,  b  and  d.  All  that 
matters  is  the  magnitudes  of  m_  and  m..  An  example  will 
illustrate  the  point. 

Consider  the  two  closed  loop  performance  functions 

pP(s)   .     1.0(B  jl) (11) 

1  (s  +  .6)2  +  (.8)2 

pp(s)  .   lOOji-rg  106)  (12) 

2  (s  +  60)2  +  (80)2 

In  both  cases  m,  »  0,6  and  nu  ■  0.8,  and  in  terms  of  T  the 
responses  are  identical. 

C(T)  -  1  +  1.118  e",6Tsin(.8T  +  6? A°   -  126. 850)    (15) 

This  response  is  shown  in  Figure  IV-3  with  three  scales 
on  the  abcissa.  The  first  scale  is  in  units  of  T.  The  second 
is  in  real  time  and  is  applicable  to  PF(s),,  where  T  =  It, 
or  tt«  T.  The  third  scale  is  also  in  real  time  and  is  appli- 
cable to  PF(s)2  where  T  ■  lOOt  or  t2=  .01T. 
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An  alternate  point  of  view  is  that  the  general  per- 
formance function  for  this  class  of  systems  has  been 
transformed  to  a  new  complex  plane  by  the  transformation 

z  -  mx  +  jm2  »  -j|p  (14) 

The  purpose  of  the  transformation  is  to  place  the  zero 
at  minus  one  on  the  z-plane.  All  possible  second  order 
systems  with  one  zero  can  be  transformed  in  this  fashion. 
If  the  variable  z  is  treated  as  a  Laplace  variable,  use 
of  the  inversion  integral  yields  an  expression  for  transient 
response,  not  in  real  time,  but  in  an  artificial  time, 
T  =  dt.  The  result  of  this  procedure  is  identical  with  the 
result  obtained  by  manipulation  of  the  expression  for  C(t). 

As  has  been  indicated,  the  transient  response  as  a 
function  of  T  is  dependent  only  on  the  values  of  m.  and  iru. 
Therefore  loci  of  constant  values  of  magnitude  of  peak 
overshoot  and  the  time  parameters  T  ,  T  ,  and  T  can  be 
drawn  on  the  z-plane. 

For  convenience  the  z-plane  may  be  regarded  as  the 
s-plane  rescaled  so  that  the  zero  is  placed  at  minus  one. 
This  places  the  complex  conjugate  poles  at  -m.  +  jm^. 
Application  of  the  Laplace  inversion  integral  yields  an 
expression  in  T  which  can  then  be  easily  converted  to  real 
time. 
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Magnitude  and  Time  of  Peak  Overshoot 

The  peak  overshoot  occurs  at  the  first  time  after  zero 
time  that  the  slope  of  the  transient  response  is  zero. 
This  requires  that 


sin(m0T  +  0)   »  0 


or 


m. 


Tr>  "  Ir^  -  h-(*   -  arctan  =^|-) 
p   m0    m0x"  1-m, 


(15) 


(16) 


Loci  of  constant  values  of  artificial  time  of  peak  over- 
shoot (T  )  are  shown  on  the  z -plane  in  Figures  IV-4  and 

rv-5. 

To  find  the  value  of  the  system  output  at  the  time  of 
peak  overshoot  we  substitute  the  expression  for  T  into 
the  expression  for  transient  response,  Equation  9* 

k2+m2  ^(l-m^+nu2  "  Sj*"-« 
1  +  _LJ: £_J ± <L-  e  d  sin(ir- 


C(Tp) 


m, 


(r-0)   (17) 


This  reduces  to 


\  m  2  +  m  2  \(l-m,)2 
C(T  )  =  1  +-l-± ' 


ra_ 


+  m^     m,  x  r 


m. 


'2      2 
e 


sin  9    (18) 


m. 


Since,  as  previously  noted,  tan  0  -  — -  ,  we  also  know  that 


m. 


sin  0  ■ 


\  mi  *  m22 


(19) 
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This  can  then  be  used  to  reduce  Equation  (18)  to 

m 

-  =i(iH*) 


2     2    m2 


C(Tp)  =  1  +  ^(1-m^)^  +  m^  e   c  (20) 

The  exponential  term  of  this  expression  is  the  amount  of 
peak  overshoot.  Magnitude  of  peak  overshoot  expressed  as 
a  percentage  of  steady  state  output  then  becomes 


M  =  100 


m 

~(*--0) 
m,  N 


\(l-m1)2  +m22  e   2     %  (21) 


Loci  of  constant  values  of  peak  overshoot  (M  )  are 

P 

shown  on  the  z -plane  in  Figures  IV-6  and  IV-7. 

It  is  of  interest  to  note  that  Equation  21  contains 

the  expression  for  peak  overshoot  of  a  pure  second  order 

ml  a 

system.  The  ratio  -—•  is  identical  to  the  ratio  ~  for  any 

given  values  of  a  and  b,  so 

e-.a/b  .  e—1/»8  (22) 

It  appears  that  the  peak  overshoot  of  a  second  order  system 
with  a  zero  is  the  peak  overshoot  associated  with  a  pure 
second  order  system  with  the  same  complex  roots,  modified 
by  a  factor  which  depends  on  the  location  of  the  zero  re- 
lative to  these  roots. 
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Rise  Time 

The  system  output  is  equal  to  its  steady  state  output 
for  the  first  time  when  the  sine  function  in  Equation  9 
is  equal  to  zero.  The  argument  of  the  sine  is  always 
negative  at  zero  time  and  increases  with  time.  Therefore 
the  first  time  after  zero  when  the  sine  function  is  zero 
is  when  the  argument  is  zero,  That  is 

m2T  +  0   -  9  =  0  (23) 


or 


I.  -  4^  (24) 


r    m2 

Loci  of  constant  values  of  artificial  rise  time  (T  )  are 
shown  plotted  on  the  z -plane  in  Figures  IV-8  and  IV- 9. 


Settling  Time 

Settling  time  is  reached  when  the  coefficient  of  the 
sine  term  in  Equation  9  decays  to  .02  or  two  per  cent  of 
the  value  of  the  steady  state  output.  This  requires  that 


1     +  m2  \(l-m1) 


,2     2 

+  nu   -m-jT 


nig 


e  A  -  .02  (25) 


A  simpler  expression  for  T  results  if  it  is  noted  that 
C(0)  =  0,  as  demanded  by  the  initial  conditions.  In  this 
case  Equation  9  becomes 


-4?- 


0  =  1  + 


-  2     2 
m1   +  m2 


\ld-i) 


2     2 


+  m 


nig 


£-  sin (0-0)        (26) 


Now  Equations  25  and  26  may  be  combined  to  yield 

"  stnOa-6)  "  -02  (27) 

From  this  the  expression  for  artificial  settling  time  is 

Ts  "  "  m"~  to  ["*02  8lnW-Q>l  (28> 

1  — 

Loci  of  constant  values  of  artificial  settling  time, 

(T_)  are  shown  on  the  z-plane  in  Figures  IV-10  and  IV-11. 

s 

The  loci  are  omitted  in  the  zone  where  peak  overshoot  is 
less  than  two  per  cent  because  settling  time  as  defined 
here  is  meaningless  in  that  area. 


Transient  Response  Parameters 
Along  the  Z -Plane  Axes 

The  values  of  the  transient  response  parameters  as 
either  m,  or  m2  approach  zero  are  of  some  interest  as 
limiting  cases.  It  is  useful  to  visualize  the  angles  0 
and  0  because  they  appear  so  often  in  the  expressions  for 
the  parameters.  They  are  shown  in  Figure  IV-12. 

First,  consider  the  parameters  as  m_  approaches  zero, 
that  is, along  the  mp  axis  or,  alternately,  on  the  s-plane 
along  the  jw  axis.  In  the  limit  we  have 
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FIGURE  IV-12 
ANGLES  9  AND  jZf  ON  THE   z -PLANE 
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T  ■  -i  (w   -  arctan  m-)  (29) 


P   m2 

2 

Up  -  100  \ 

1  +  m22  ^ 

m     m  -   arctan  m  +  %/2 

r          nu 
2 

T  is  infinite, 
s 

(30) 
(3D 


As  ra2  approaches  zero  there  are  two  distinctly  dif- 
ferent possibilities.  For  nu  less  than  one,  in  the  limit 
as  m2  approaches  zero,  T  ,  T  ,  and  T  become  infinite  and 
M  becomes  zero. 

On  the  other  hand  if  nu  is  greater  than  one,  in  the 
limit  as  iru  approaches  zero,  we  have 


TP  "  "  l-*i      ^  (52) 

.  +  1-m,  oi 
Mn  -  100  (1-m  )  e     L     /o  (33) 

P  J. 

T  is  infinite, 
s 

These  expressions  are  useful  in  determining  the  in- 
tercepts of  a  particular  locus. 
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Construction  of  the  Loci 

Each  of  the  families  of  loci  were  constructed  by  cal- 
culating and  plotting  the  values  of  the  parameters  as 
functions  of  m2  for  a  number  of  fixed  values  of  m..  The 
result  was  a  series  of  cross-sections  of  the  surfaces  de- 
fined by  the  parameters  over  the  z-plane.  The  coordinates 
m,  and  m^  were  then  read  from  these  cross -sect ion  curves 
for  each  locus  desired,  and  these  points  were  plotted  on 
the  z-plane.  This  method  was  required  because  no  method 
could  be  found  to  determine  the  equations  of  the  loci  on 
the  z-plane. 

Use  of  the  Curves 

The  curves  are  used  to  determine  the  transient  response 
parameters  in  the  same  manner  as  outlined  in  Chapter  III 
in  connection  with  pure  second  order  systems.  Once  the 
closed  loop  performance  function  is  known  in  factored 
form,  that  is,  the  quantities  a,  b,  and  d  are  determined, 
m.  and  nu  must  be  calculated.  This  is  relatively  simple 
since 

■i-iSr   md    m2"A 

The  point  -m,  +  Jnu  is  then  located  on  each  family  of 
loci  and  values  of  T  ,  Tr,  Ts,  and  M  are  read  off  directly. 
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The  three  time  parameters  can  then  be  converted  to  real 
time  by  dividing  each  by  the  value  of  d,  i.e., 


T  T  T 

t     =  t§t-,  t     =  ttt>         and     t„  =  -rlr 

p        |d|  '  r        |d|  '  s        |d| 


With  these  three  times  and  the  magnitude  of  peak  over- 
shoot the  transient  response  may  be  rapidly  sketched  if  it 
is  remembered  that: 


a.  The  slope  at  t  -  t  is  zero. 

P 

b.  The  transient  oscillating 
frequency  is  -^   cycles  per 
second. 

c .  The  sinusoid  is  zero  at 

t  =  t  ,  and  is  increasing 
with  time. 

d.  The  slope  at  zero  time  is 

2  2 
a  +b 

T^  ,  which  can  be  easily 

shown  by  use  of  the  initial 

value  theorem. 


Observations 


The  loci  of  constant  values  of  peak  overshoot  become 
straight  lines  as  they  approach  the  origin.  If  iru  and  m 
are  much  less  than  one,  the  expression  for  peak  overshoot- 
degenerates  to  m. 


% 

m0 
M  -  100  e   2 
P 
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As  has  been  shown,  this  is  identical  to  the  expression  for 
peak  overshoot  for  pure  second  order  systems.  These  loci 
are  actually  very  nearly  straight  lines  if  m,  and  m^  are 
less  than  about  0.6,  indicating  that  a  zero  will  not  in- 
validate a  pure  second  order  approximation  unless  its 
magnitude  is  less  than  about  twice  the  system  undamped 
natural  frequency.  For  very  low  values  of  both  m,  and  m. 
the  loci  of  constant  values  of  parameters  have  been  omitted 
because  the  pure  second  order  curves  are  very  nearly  exact 
in  this  region. 

A  comparison  of  the  loci  of  constant  values  of  peak 
overshoot  shown  in  Figure  IV-7  with  the  same  loci  for  pure 
second  order  systems  as  shown  in  Figure  III-3  will  illustrate 
the  previous  observation.  The  comparison  will  also  amplify 
the  phrase  "relatively  more  oscillatory  than  would  be  the 
case  without  the  zero".  For  fixed  complex  root  locations 
the  addition  of  a  zero  can  cause  a  large  increase  in  the 
magnitude  of  peak  overshoot.  In  fact,  peak  overshoot  of 
greater  than  one  hundred  percent  is  a  distinct  possibility 
for  such  systems .  It  should  be  noted  that  the  addition  of 
a  zero  does  not  effect  the  transient  oscillating  frequency. 

The  initial  slope  of  the  transient  response  is  not  zero 
as  was  the  case  with  the  pure  second  order  systems .  The 
initial  conditions  are  met,  however,  because  there  is  a 
velocity  discontinuity  at  zero  time  due  to  the  combination 
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of  the  input  function  and  the  zero.  As  was  mentioned, 
the  zero  can  arise  from  the  addition  of  derivative 
control  in  the  forv/ard  path.  Differentiation  of  the 
unit  step  through  this  channel  yields  a  unit  impulse, 
and  it  is  this  impulse  that  accounts  for  the  initial 
slope. 
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Chapter  V 
THE  PURE  THIRD  ORDER  SYSTEM 

The  addition  of  a  real  root  to  the  pure  second  order 
system  results  in  a  new  system  which  we  define  as  a  pure 
third  order  system.  The  closed  loop  performance  function 
of  the  pure  third  order  system  then  has  one  pair  of  com- 
plex conjugate  poles,  one  real  pole  and  no  zeros. 

Much  less  is  known  about  the  pure  third  order  system 
than  is  the  case  with  the  pure  second  order  system.  This 
is  due  largely  to  the  introduction  of  a  third  variable  in 
the  equations  due  to  the  third  root.  The  mathematical 
expression  for  the  transient  response  of  the  pure  third 
order  system  contains  a  decaying  exponential  term  in  ad- 
dition to  the  constant  and  damped  sinusoidal  terms  of  the 
second  order  system.  The  actual  expression  for  the 
transient  response  is  not  difficult  to  obtain,  nor  to  plot, 
but  rather  the  actual  plotting  is  laborious .  The  main 
difficulty  in  establishing  standard  sets  of  curves  for  the 
third  order  system,  as  is  often  done  with  the  second  order 
system,  is  in  portraying  these  curves  as  functions  of 
three  variables. 

This  does  not  imply  that  nothing  has  been  done  in 
this  area.  To  the  contrary,  much  effort  has  gone  into  the 
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study  of  third  order  systems.  Unfortunately,  the  results 
of  this  effort  has  often  produced  new  ways  of  factoring 
the  system  characteristic  equations  without  producing  any 
simple  method  of  evaluating  the  transient  response.  One 
such  method  was  developed  by  Y.  J.  Liu  (5),  and  is  based 
on  nondimensionalizing  the  characteristic  equation.  By 
using  this  method,  stability  and  the  roots  of  the  char- 
acteristic equation  can  be  determined  from  a  set  of  non- 
dimensional  i  zed  charts. 

Another  method  of  importance  in  any  discussion  of 
third  order  systems  is  Mitrovic ■ s  method  (6 ) .  In  this 
method,  by  use  of  a  selected  transformation,  it  is 
possible  to  compute  a  set  of  curves  which  inherently  con- 
tain the  solution  to  all  third  order  equations  and  can  be 
used  to  determine  stability.  These  curves  have  been  com- 
puted and  are  known  as  third  order  charts.  They  may  be 
used  directly  in  compensation  by  cancellation  provided 
suitable  root  locations  are  known.  However,  this  still 
does  not  solve  the  problem  of  where  to  locate  the  roots 
and  it  does  not  yield  transient  response  directly. 

Transient  Response 

A  plot  of  the  roots  of  the  pure  third  order  system  is 
shown  in  Figure  V-l.  The  closed  loop  performance  function 
for  the  unity  feedback  third  order  system  may  be  written  as 
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FIGURE  V-l 
POLES  OP  THE  PURE  THIRD 
ORDER  SYSTEM  ON  THE  s -PLANE 
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PF 


(3)  ~  Rlli  ■  (s4c)(s+a+jbj(s+a-jb)  (la) 


This  may  be  simplified  as  follows 

PF<S)  -  ;    (?%/b2ii — rr  (lb) 

(s-K>)  [(s+ar  +  b^J 

If  a  unit  step  input  is  applied  to  a  system  having 
the  performance  function  shown  in  Equation  1,  the  complex 
frequency  domain  output  becomes 

„(,),    <»2+b2)% _-  (2) 

s(s+c)  |_(s+a)  +  b  J 

Th®  inverse  of  this  transform  gives  the  time  domain  out- 
put of  the  Bystem,  which  is 


C(t)  -  1  - 


(a2  +  b2)   -ct 


fta-c)2  +b2J 

dZTEL  e-at  8in(bt_0.^)     (3) 


b  ]|(c-a)2  +  b2 


where 


0  ■  arctan  — - 
-a 


fi  -  arctan  — - 


Now  if  we  redefine  the  root  locations  in  terms  of  the 
location  of  the  real  root  as  shown  in  Figure  V-2,  the  real 
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ILLUSTRATION  OP  QUANTITIES  n.,  AND  n 
PURE  THIRD  ORDER  SYSTEM 
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and  imaginary  parts  of  the  complex  roots  become 


ft  •  H_  |C 


n2|c 


The  quantities  n,  and  n?  are  merely  rescaled  values  of  the 
real  and  imaginary  parts  of  the  complex  roots. 

Making  these  substitutions  for  a  and  b  in  Equation  5, 
the  expression  for  transient  response  becomes 

2     2 

C(t)  =  1 1  g  2  g  e"ct 

(l-n^*  +  n2* 


\  n_   +  n0        -n.  ct 

e  -1  sin(noct-9-0)      (4) 


n 


where 


2  ^(l-n^2  +  n22 


9  =  arc  tan  i^/-bl 

n2 
0   =  arc  tan  -n~- 
^         1-n, 


This  substitution  has  partially  nondimensionalized 
the  transient  response  of  systems  with  the  same  relative 
root  configuration.  The  transient  response  can  now  be 
completely  nondimensionalized  by  making  the  additional  sub- 
stitution T  m    |c|te  This  is  equivalent  to  defining  a  new 
artificial  time  for  this  class  of  systems  as  was  done  in 
Chapters  III  and  IV.  The  transient  response  now  becomes 
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(n  2  +  n  2)   _T 
C(T)  =  1       .&     g  e 


(l-n^*  +  n2' 


(5) 


\  n,  +  n, 


-n..T 

.     x     sin(noT-Q-0) 


n2  \j  (1-n^2  +  n. 


It  is  of  interest  to  look  at  a  different  derivation 
of  Equation  5.  As  a  first  step  let  us  write  the  closed 
loop  performance  function  in  terms  of  the  root  locations 
shown  in  Figure  V-2, 


PF(s) 


(nx  +  n2  )c^ 
(s+c)  [(s-ttyO  +  (n2c)cJ 


Dividing  numerator  and  denominator  by  c  gives 


PF(s 


2     2 
nl  +n2 


(S  +  1)  [(2  +  nx)2  +  n. 


Now  if  we  make  the  linear  transformation 


Z  -  «!  +  Jn2  "  "lei 


(6a) 


(6b) 


the  performance  function  becomes 


PF(z) 


2     2 
nl  +n2 


(z+l)  [(z-ii^)  +  n2 


(7) 


This  has  effectively  transformed  the  system  onto  a 
new  complex  plane  in  which  the  real  root  is  located  at 
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minus  one  and  the  complex  roots  are  at  -n_  +  jn  .  Now 
upon  application  of  a  unit  step  the  complex  frequency 
domain  output  becomes 

r,     2        2 

n   +  n0 

C(z) 1     r  2 2 «  (8) 

z(z+l)   (z-fn, )  +  np  J 

The  inverse  transform  of  this  gives  the  transient  response 


2     2 
nn   +  n0      « 

C(T)  -  1 s  2 ^  e  i 

(1-n..  )  +  n0c 

1      2  (9) 


+ 


i   2      2 

\n,   +  n0        -n,T 

e  x   sin(noT-9-0) 


n2  ^(1-n^2  +  n22 


In  this  case  T  is  the  time  associated  with  the  inverse 
from  the  z-plane.  Note  that  Equation  9  is  identical  to 
Equation  5. 

Time  and  Magnitude  of  Peak  Overshoot 

As  in  the  previous  cases,  the  slope  of  the  transient 
response  is  zero  at  the  time  of  peak  overshoot.  That  is 

&rC(tp)-0  (10) 

The  Laplace  transform  of  this  for  zero  initial  conditions 
is 
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roc(*) 


-  S  C(s) 


(11) 


Multiplying  both  sides  of  Equation  2  by  s  gives 


S  C(s) 


(a2  +  b2)c 


(s+c)  [(s+a)2  +  b2] 


(12) 


Combining  Equations  11  and  12,  and  obtaining  the  inverse 
transform  gives  the  slope  of  the  transient  response  as 


a  „(tJ.  (a2  y2>«  >"otP 

(c-a)2  +  b2 


3E  C* 


+     (a2  ±  b2) 


b  \|(c-a)w 


-  e   P  sin(bt  -0) 
2  p 


b_ 
c-a 


where  0   is  defined  as  before    0  ■  arctan  — 
Equation  15  can  be  reduced  to 


(13) 


-at 


(a2+  b2)[b  e"°  P+  ^(c-a)2+b2  e"&  Psin(btp  -  0)]  m  Q  (l4) 

b  [(c-a)2  +  bS 


For  practical  systems  the  left  side  of  Equation  14  can 

equal  zero  only  when  the  term  in  brackets  equals  zero. 

Therefore,  the  condition  to  be  satisfied  at  time  of  peak 

overshoot  is 

-ct     I    o    ^  -at 
be   p  +  \|(c-a)2  +  b2  e   p  sin(btp  -  0)   -  0       (15) 
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Equation  15  can  now  be  nondimensionalized  by  making  the 
same  substitutions  as  used  in  the  transient  response,  i.e., 


a 
b 
T 


n-Jc 

n2|c 

Iclt 


When  these  substitutions  are  made  in  Equation  15  the  con- 
dition to  be  satisfied  at  time  of  peak  overshoot  becomes 


-n, 


-(l-n-^T 


3in(n0T  -  0) 


(16) 


(1-n^  +n2' 


n2 

where    0  ■  arc  tan  -r-~ 

l-n 


Prom  Equation  5  we  see  that  the  magnitude  of  peak 
overshoot  (expressed  as  a  percentage)  can  be  written  as 


M_ 


100 


~  2     2 
"l   +n2 


n. 


(l-i^)1 


+  n 


-n.T 
—  e  x   p  sin(n0T  -9-0) 

2  *   P 


2     2 
n,   +  n^ 


2 

-T 


(1-^)' 


+  n, 


7. 


(17) 


For  any  combination  of  n.  and  n_,  Equation  16  can  be 
solved  by  trial  and  error  techniques  to  determine  the 
values  of  T  for  the  first,  second  and  each  succeeding  maxi- 
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mum  of  the  transient  response.  With  a  small  amount  of 

practice  these  values  of  T  can  be  determined  with  one  to 

three  iterations  per  value.  The  values  of  T  determined 

P 

from  Equation  16  can  then  be  inserted  in  Equation  17  to 
determine  which  maximum  gives  the  peak  overshoot,  and 
the  amount  of  that  overshoot.  From  the  foregoing  we  see 
that  the  calculations  for  time  and  magnitude  of  peak  over- 
shoot are  inseparable.  Each  must  be  known  before  the 
other  can  be  determined  correctly.  This  is  not  a  severe 
handicap,  however,  as  the  peak  overshoot  occurs  on  the 
first  maximum  for  practically  all  combinations  of  n_  and 
n2  which  include  values  of  n_  less  than  two. 

Loci  of  constant  values  of  magnitude  of  peak  overshoot 
(M  )  are  shown  on  the  z -plane  in  Figure  V-3.  Another  set 
of  loci  of  constant  values  of  peak  overshoot  with  a  com- 
pressed ordinate  scale  to  permit  extension  of  the  plot  to 
include  larger  values  of  n2  is  shown  in  Figure  V-4.   Although 
Figure  V-4  does  not  represent  a  complex  plane  due  to  the 
difference  in  scales  used  on  the  axes,  it  may  be  considered 
as  one  for  the  purpose  of  reading  values  of  peak  overshoot 
associated  with  values  of  n,  and  n  . 

Loci  of  constant  values  of  artificial  time  of  peak 
overshoot  (T  )  are  shown  on  the  z-plane  in  Figure  V-5. 
An  interesting  point  in  connection  with  T  is  that  as  n2 
increases  the  peak  overshoot  occurs  on  the  second  or  sub- 
sequent maximum  of  the  transient  response  curve,  and  the 
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time  of  peak  overshoot  is  no  longer  unique.  The  zone  in 
which  this  situation  exists  is  shown  in  Figure  V-6  which 
is  an  extension  of  Figure  V-5. 

Rise  Time 

Rise  time  is  the  first  time  at  which  the  output  of  a 
system  equals  its  steady  state  value.  This  is  equivalent 
to  saying  that  rise  time  is  the  first  time  for  which  the 
overshoot  is  zero.  From  Equation  5  we  see  that  the  con- 
dition for  zero  overshoot  is 


2     2 


2     2 


nn   +  n9      ™    \n   +  n0  -nnT 

-= *-= w  e   =  '  x     g       e  x 


.2  .  _  S 


(l-i^)   +n2       ^  ^(1.ni)2  +  ^ 
This  can  be  rearranged  as  follows 


2 


e  1  sin(noT-9-0)    (18) 


ln(n2T=9-0)  -  n, 


nn  2  +  n^2      - (1-n, )T 

e     X  (19) 


2     2 


(1-n^)  +  n2 


Obviously,  this  too  is  a  transcendental  equation,  but 
just  as  in  the  case  of  time  of  peak  overshoot,  it  can  be 
solved  by  trial  and  error  methods  for  any  selected  com- 
bination of  n,  and  n2. 

Loci  of  constant  values  of  artificial  rise  time  (T  ) 

r 

are  shown  on  the  z -plane  in  Figure  V-7. 
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Settling  Time 

As  normally  defined,  settling  time  is  the  last  time  at 
which  the  output  error  exceeds  some  fixed  percentage  of 
the  steady  state  output.  This  error  is  often  specified  as 
two  percent.  Using  this  definition,  settling  time  is  ex- 
tremely difficult  to  determine  for  a  third  order  system  as 

the  condition  that  must  exist  at  T  is 

s 


2      2  2      2 

n,   +  n0  -n-T  n.c  +  n0      „, 

e  x   sin(noT-0-j2f)   — ^ — r^ ^  e"1  -  +  .02 

+  nrt 


>  |  (l-nx)2  +  n22  (lHV 

(20) 


To  determine  settling  time  exactly  would  require  several 
trial  and  error  solutions  of  Equation  20  or  a  plot  of 
Equation  5  for  each  combination  of  n.  and  n_.  The  labor 
involved  in  either  method  is  prohibitive. 

The  settling  time  for  the  third  order  system  as  defined 
here  is  the  time  at  which  the  sum  of  (1  -  exponential  term) 
plus  the  height  of  the  envelope  of  the  sinusoidal  term  is 
equal  to  two  percent  of  the  steady  state  output.  The 
equation  to  be  satisfied  for  this  condition  is 


2     2  2     2 

n.   +  n0       -niTa     n,   +  n0     ~T« 

e  ls+ 1  g  2  g  e  s  -  .02    (21) 


n2 


(l-^)2  +n22         (W  +  n2 


Loci  of  constant  values  of  artificial  settling  time 

(T  )  are  shown  on  the  z-plane  in  Figure  V-8.  Settling  time 

s 
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for  systems  with  less  than  two  percent  overshoot  is  not 
shown  because  settling  time  as  defined  here  is  meaningless 
for  such  systems. 

Construction  of  the  Loci 

In  constructing  the  loci  of  constant  parameter  values 

it  was  necessary  to  calculate  the  values  of  each  parameter 

for  a  large  number  of  combinations  of  n-  and  np.  This  was 

by  calculating  the  parameter  values  as  a  function  of 


np  for  a  number  of  fixed  values  of  n_.  All  of  the  calcu- 
lations were  performed  by  trial  and  error  as  previously 
indicated. 

The  values  of  the  parameters  were  then  plotted  as  a 
function  of  n,  or  n_  with  the  other  variable  n_  or  np  as 
a  parameter  of  the  curves.  This  gave  a  series  of  cross- 
sections  of  the  surfaces  defined  by  the  transient  response 
parameters  over  the  z-plane.  The  coordinates  n,  and  n~ 
were  then  read  from  these  cross-section  curves  for  each 
locus  desired.  These  points  were  plotted  on  the  z-plane 
and  the  constant  value  loci  drawn. 

Observations 

The  loci  of  constant  values  of  peak  overshoot  become 
straight  lines  as  they  approach  the  origin  of  the  complex 
plane.  This  is  consistent  with  the  usual  second  order 
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approximation  for  the  case  in  which  the  value  of  the  real 
pole  is  several  times  the  value  of  the  undamped  natural 
frequency  of  the  system.  In  Figure  V-9  the  validity  of 
this  approximation  is  demonstrated  by  a  comparison  of  the 
peak  overshoot  of  the  pure  second  order  system  with  that 
of  the  pure  third  order  system  for  small  values  of  n-  and 

v 

Another  area  of  interest  on  the  plot  of  overshoot  loci 
is  the  area  of  zero  overshoot.  The  zero  overshoot  locus 
represents  the  boundary  where  the  magnitude  of  peak  over- 
shoot becomes  too  small  to  calculate.  The  boundary  as 
shown  represents  a  value  of  the  order  of  e~   .  For  values 
of  n,  greater  than  about  0.94  the  peak  overshoot  is  too 
small  to  calculate,  regardless  of  the  value  of  n2.  In  this 
case  the  response  is  essentially  that  of  a  first  order 
system  with  a  decaying  ripple  superimposed  on  the  rising 
exponential  curve. 

As  noted  previously,  there  exists  a  region  in  which 
the  time  of  peak  overshoot  is  not  necessarily  unique. 
This  area  is  shown  in  Figure  V-6.  The  zone  is  bounded  by 
a  line  through  the  values  of  i^  and  n^  for  which  the 
height  of  the  first  and  second  maximums  of  the  transient 
response  curve  are  identical.  Obviously  this  includes 
the  entire  imaginary  axis. 

For  small  values  of  n2  the  time  of  peak  overshoot  is 
relatively  independent  of  the  value  of  n,.  In  fact,  the 
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only  area  where  time  of  peak  overshoot  becomes  sensitive 
to  the  value  of  n,  is  where  both  n..  and  n?  are  large. 

For  small  values  of  n2  the  rise  time  is  also  relatively 
independent  of  the  value  of  n,,  although  it  is  more  dependent 
on  this  value  than  the  time  of  peak  overshoot  is.  For  large 
values  of  both  n1  and  n-  the  rise  time  becomes  more  depend- 
ent on  the  value  of  n1  than  on  the  value  of  n_. 

For  small  values  of  n.  the  settling  time  is  essentially 
independent  of  the  value  of  n.,  and  settling  time  is  never 
as  dependent  on  the  value  of  n^  as  it  is  on  the  value  of  n_. 

To  avoid  heavy  concentrations  of  lines  the  loci  are  not 
plotted  near  the  origin  of  the  complex  plane.  This  does 
not  detract  from  their  usefulness  as  the  second  order 
approximation  is  valid  in  this  area. 

Use  of  the  Curves 

The  curves  for  the  pure  third  order  system  are  used  in 
exactly  the  same  manner  as  the  curves  of  Chapters  III  and 
IV.  Since  their  use  has  already  been  explained  there, 
the  explanation  will  not  be  repeated  here.  Instead  an 
actual  example  is  presented. 

Assume  that  a  system  has  the  following  closed  loop 
performance  function 

2500 


K    '   "(s-f5)(s+2+j4)  (B+2-J4) 


•77- 


Dividing  the  real  and  imaginary  parts  of  the  complex  poles 
by  the  value  of  the  real  pole  gives: 

2    u 
•l  "  f  "  -4 

n2  -  }  -  .8 

Entering  Figures  V-3,  V-5,  V-7,  and  V-8  with  these  values 
of  n,  and  n2  we  can  read  the  following  parameter  values: 

Mp  -  1% 

Tp  -  5.1 

Tr  -  3.8 

T  =  10,2 
s 

Dividing  the  three  artificial  times  by  the  value  of  the 
real  pole  gives  the  real  times: 

t  -  %^     -  1.02  sec. 
p    5 

t  -  2j£    .  ,76  sec. 
r    5 

t  =  ^2^    -  2.04  sec. 
s    5 

From  the  performance  function  we  see  that  the  transient 
oscillating  frequency  is  four  radians  per  second,  or  the 

IT 

period  is  p  seconds.  The  three  values  of  time,  the  period 
of  oscillations,  magnitude  of  peak  overshoot  and  the  fact 
that  the  slope  is  zero  at  zero  time  and  at  time  of  peak 
overshoot  constitute  sufficient  information  from  which  to 
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sketch  the  transient  response.  The  transient  response  is 
sketched  in  Figure  V-10. 

Now  if  we  apply  a  unit  step  input  to  the  system  and 
obtain  the  transient  response  by  performing  the  inverse 
transformation,  we  get 

C(t)  =  1  -  8e"5t  +  l.ll6e"2tsin(229t  -  169.7) 

This  output  has  been  calculated  and  is  plotted  as  a  function 
of  time  in  Figure  V-ll.  From  an  expanded  scale  plot  of  this 
output,  transient  response  parameters  have  been  determined. 
These  are  compared  with  the  predicted  values  below. 

Predicted  Value  from 

Parameter  Value  Transient  Response 

t  1.02  sec.  1.025  sec. 

P 

t  .76  sec.  .76  sec. 

t  2.04  sec.  2.00  sec. 

s 

Mp  13  #  12.7  % 
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Chapter  VI 
EXTENSION  OF  THE  TECHNIQUE 

The  curves  previously  developed  can  be  used  to  approxi- 
mate the  transient  response  of  third  order  systems  with 
one  or  two  zeros.  This  extension  is  made  possible  by  a 
simple  algebraic  manipulation  which  breaks  a  system  down 
into  two  or  three  parallel  systems.  The  transient  response 
of  each  of  the  parallel  systems  can  be  approximated  using 
the  curves  of  Chapters  III,  IV  and  V.  These  approximate 
transient  responses  can  then  be  combined  to  give  an  approxi- 
mation of  the  transient  response  of  the  original  system. 

The  scheme  can  be  explained  best  by  an  illustration. 
Consider  the  closed  loop  performance  function 

PP(B)  .  c(a2+b2) is^ (1) 

*      (s+c)  [(s+a)2  +  b2J 

This  can  be  rewritten  as 

PP(s)  -  e(a2+b2)v    (f+c+f)  (2) 

d       (s+c)  [(s+a)2  +  b2J 

where  f  is  defined  by 

f  -  (d  -  c)  (5) 
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We  can  then  write  the  performance  function  as 


PF(s) 


c 
d" 


o2     H2 

a  +  t> 


(s+a)2  +  b2 


♦§ 


(a2  +  b2)c 


(s+c)[(s+a)2  +  b2J 


(4) 


Applying  a  unit  step  input  and  taking  the  inverse  transform 
gives  the  transient  response  as 


c(t)  -# 


Pure  second 
order  response. 


+  £ 


Pure  third 
order  response 


To  sketch  this  transient  response  the  pure  second  order 
response  and  the  pure  third  order  response  must  both  be 
sketched  first  by  the  methods  previously  described.  Both 
must  be  sketched  as  functions  of  real  time.   At  any  instant 

the  total  response  is  §■  multiplied  by  the  value  of  the  pure 

f 
second  order  response  at  that  time  plus  -r  multiplied  by  the 

value  of  the  pure  third  order  response  at  the  same  time. 

The  two  responses  can  be  combined  most  rapidly  by  graphical 

methods . 

As  an  example,  consider  the  performance  function 


PF(s) 


26(s+2) 


(s-*4)[(s+2)2  +  52] 
In  this  case     f  -    (d  -  c)  ■    (2  -  4)  -  -2.     Therefore 


(5) 


PF(s)  -  § 


13 

~7S 


(s+2r  +  5* 


2 
2 


52 


(s44)[(s+2)2  +32] 


(6 


C(t) 


Pure  second 
order  responsej 


-1 
-83- 


Pure  third 
order  response 


For  the  pure  second  order  response,  a  =  2  and  b  *  3 
From  the  curves  in  Chapter  III  values  of  parameters  are 


TP 

■ 

2, 

1 

Tr 

- 

1. 

,47 

T» 

- 

4, 

1 

M 

. 

12 

>.8# 

p 

Dividing  the  artificial  times  by  the  value  of  a  yields 


*p 

= 

1 

.05 

sec. 

*r 

s 

.735 

sec. 

*a 

- 

2 

.05 

sec. 

For  the  third  order  response: 


nl  *  T\   *  '-> 


n2  *  4  '    *75 


From  the  curves  in  Chapter  V  values  of  parameters  are 


TP 

3B 

5. 

,6 

Tr 

m 

4, 

.3 

T» 

S= 

8, 

.5 

MP 

" 

6, 

.* 
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Dividing  these  artificial  times  by  the  value  of  c  yields 


t  =  1.4    sec. 

P 

t  »  1.075  sec. 


t  -  2.125  sec. 
s 


Note  that  in  both  cases  the  transient  oscillating  frequency 

■5 
is  4-  cycles  per  second. 

The  process  of  combining  the  two  responses  is  accomplished 
graphically.  This  is  illustrated  in  Figure  VT-1  where  the 
two  transient  response  curves  and  their  sum  have  been  sketched 

The  same  method  can  be  used  to  approximate  the  response 
of  third  order  systems  with  two  zeros.  This  merely  involves 
more  labor  and  of  course  the  accuracy  declines  with  the  ad- 
dition of  zeros. 

Consider  the  following  closed  loop  performance  function 

PPfB)  -  o(ag+b2)   (s^Ha+d  ) 

ala2     (s+c)[(s+a)2  +  b2] 

This  can  be  rewritten  as 


dld2     (8+e)  [(s+a)2  +  b2] 


(8) 


where 


f  1  "  *1.  "  ° 

f 2   »  d2  -  c 


-85- 


Multiplying  the  terms  in  the  numerator  yields 

PF(3),cj4^,(^>g  +  (vv(-)+v2      {9) 

dld2       (s-kj)  [(s+a)2  +  b2] 


This  can  be  broken  into  three  terms  as  follows : 


PWO    °2   a2+b2     s^       ,  (fl+f2>c    a2+b2 

(10) 

,  flf2     c(a2  +  b2) 
dld2   (s-fc)[(s+a)2  +  b2] 


This  can  be  sketched  in  the  same  manner  as  before  except 
that  the  transient  response  curve  will  be  the  sum  of 
three  curves  instead  of  two. 

No  satisfactory  method  of  approximating  the  transient 
response  of  more  complicated  systems  has  been  found.  The 
only  alternatives  seem  to  be: 

a.  Approximate  higher  order  systems 
with  lower  order  systems  by 
ignoring  the  poles  and  zeros  that 
are  far  out  in  the  left  hand  plane, 
or 

b.  Expand  the  performance  function  by 
the  method  of  partial  fractions  into 
fragments  of  manageable  size,  utilizing 
the  techniques  developed  in  this  thesis 
whenever  possible. 
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Chapter  VII 
CONCLUSIONS  AND  RECOMMENDATIONS 

The  servomechanisms  investigated  in  this  thesis  have 
been  divided  into  classes  according  to  the  number  and  type 
of  poles  and  zeros  of  their  closed  loop  performance 
functions.  This  was  done  to  facilitate  the  investigation 
of  the  transient  responses  of  the  different  classes  of 
systems.  Systems  investigated  in  this  thesis  are  those 
whose  closed  loop  performaance  functions  contain  the 
following  pole-zero  configurations: 

a.  Two  complex  poles  only. 

(Pure  second  order  systems) 

b .  Two  complex  poles  and  one 
zero.   (Second  order  system 
with  one  zero) 

c.  Two  complex  poles  and  one 
real  pole.  (Pure  third 
order  system) 

It  has  been  shown  that  for  each  of  these  classes  of 
systems  the  transient  response  of  systems  with  the  same 
relative  pole-zero  configuration  can  be  nondimensionalized 
The  transient  response  can  be  expressed  as  a  function  of 
an  artificial  time  and  two  complex  coordinates.   For  the 
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pure  second  order  system  the  artificial  time  is  the  magni- 
tude of  the  real  part  of  the  complex  poles  multiplied  by 
real  time  and  the  complex  coordinates  are  the  real  and 
imaginary  parts  of  the  complex  poles.   For  systems  with 
three  singularities  the  artificial  time  is  the  magnitude 
of  the  real  pole  or  zero  multiplied  by  real  time  and  the 
complex  coordinates  are  the  real  and  imaginary  parts  of 
the  complex  poles  divided  by  the  magnitude  of  the  real 
pole  or  zero. 

With  the  transient  response  expressed  in  this  manner, 
magnitude  and  time  parameters  of  the  response  can  be 
expressed  in  terms  of  percentage  of  steady  state  output 
and  artificial  time  respectively.  Loci  of  constant  values 
of  four  transient  response  parameters  for  each  of  the 
three  classes  of  systems  listed  above  have  been  mapped  onto 
a  complex  plane.  These  parameters  are  magnitude  and  arti- 
ficial time  of  peak  overshoot,  artificial  rise  time  and 
artificial  settling  time.  For  the  pure  second  order  system 
the  loci  are  mapped  directly  onto  the  s -plane,  while  for 
the  systems  with  three  singularities  the  loci  are  mapped 
onto  an  auxiliary  complex  plane. 

The  approximate  value  of  these  four  parameters  for 
practically  any  specific  system  of  the  three  classes  dis- 
cussed can  be  determined  from  these  families  of  loci  of 
constant  parameter  values.  A  rather  accurate  transient 
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response  curve  can  then  be  sketched  using  these  four  para- 
meters, transient  oscillating  frequency,  initial  slope  and 
slope  at  time  of  peak  overshoot.  This  has  been  demonstrated 
by  example. 

Although  the  discussion  in  this  thesis  has  been  limited 
to  systems  which  exhibit  unity  feedback  characteristics, 
this  is  not  a  necessary  restriction.  The  closed  loop  per- 
formance function  of  a  non-unity  feedback  system  can  be 
written  as  a  constant  multiplied  by  the  performance  function 
of  a  unity  feedback  system.  Regardless  of  the  value  of  this 
constant,  it  will  appear  in  the  transient  response  as  a 
multiplier  of  all  terms.  Therefore,  the  value  of  the  con- 
stant will  not  effect  the  value  of  any  of  the  transient 
response  parameters  determined  in  this  thesis.  The  transient 
response  of  such  a  system  will  merely  have  a  steady  state 
value  different  from  unity. 

Although  loci  of  only  four  transient  response  parameter® 
have  been  mapped  onto  complex  planes  in  this  thesis,  loci 
of  other  parameters  could  be  mapped  in  a  similar  manner. 
Other  such  parameters  which  would  appear  to  be  useful  in 
improving  the  resulting  transient  response  sketches  include? 

a.  Delay  Time  (Time  at  which  the  output 
first  equals  one-half  the  steady 
state  value). 
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b.  Time  and  Magnitude  of  First 
Undershoot . 

Inclusion  of  additional  singularities  in  the  closed 
loop  performance  funtion  of  a  system  presents  additional 
problems.   Although  the  technique  developed  in  this  thesis 
can  be  applied  to  such  systems,  one  set  of  loci  is  not 
sufficient  for  all  systems  of  a  class.  Consider  a  system 
with  two  complex  and  two  real  singularities.  Application 
of  the  mapping  technique  would  require  fixing  one  of  the 
real  variables  as  a  function  of  the  other  real  variable 
and  construction  of  a  set  of  loci  for  this  fixed  relative 
configuration  of  the  two  variables.  Additional  families 
of  loci  would  be  required  for  each  fixed  configuration 
of  these  two  real  variables  for  which  knowledge  of  the 
parameters  was  desired. 

One  method  of  avoiding  the  additional  mapping  required 
for  third  order  systems  with  one  or  two  zeros  has  been 
presented.  This  involves  separation  of  such  systems  into 
parallel  components  whose  individual  responses  can  be 
approximated  using  the  parameter  loci  presented  in  this 
thesis.  Once  the  component  responses  have  been  approxi- 
mated they  may  be  combined  to  give  a  resultant  approxi- 
mation of  the  system  transient  response.  Using  this 
scheme  the  accuracy  of  the  resultant  approximation  is 
considerably  less  than  the  accuracy  of  the  component 
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approximations.   However,  careful  construction  of  the  com- 
ponent approximations  should  permit  the  construction  of  a 
useful  resultant  approximation,  at  least  for  the  third 
order  system  with  one  zero. 

The  technique  developed  in  this  thesis  is  not  restricted 
to  the  mapping  of  transient  response  parameters.   Frequency 
response  parameters  could  also  be  mapped  onto  a  complex 
plane  by  a  similar  technique. 

The  authors  acknowledge  that  they  have  not  exhausted 
the  possibilities  of  the  technique  presented  in  this 
thesis.   In  fact,  they  believe  that  there  are  a  vast  number 
of  possibilities  yet  to  be  investigated.  However,  due  to 
the  labor  involved  in  the  mapping  technique,  further  investi- 
gations should  be  conducted  with  the  aid  of  a  digital  com- 
puter.  Specific  areas  recommended  for  investigation  using 
a  computer  ares 

a.  The  three  classes  shown  in  this 
thesis  to  extend  the  loci  plots 
to  include  larger  values  of  the 
complex  coordinates. 

b,  The  third  order  system  with  one 
zero  to  include  construction  of 
sets  of  parameter  loci  for  each 
of  a  number  of  fixed  relative 
configurations  of  the  real  pole 
and  zero. 
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c.  Construction  of  loci  of  ad- 
ditional parameters  such  as 
delay  time  and  time  and 
magnitude  of  the  first  under- 
shoot for  the  systems  listed 
in  a  and  b. 

d.  Construction  of  loci  of 
constant  values  of  frequency 
response  parameters. 

e.  The  possibility  of  mapping 
statistical  quantities  such 
as  integral  error  squared 
for  various  inputs. 
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